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Abstract 



We analyze the dynamics of the Mixmaster Universe on the base of a standard 
Arnowitt-Deser-Misner Hamiltonian approach showing how its asymptotic evolu- 
tion to the cosmological singularity is isomorphic to a billiard on the Lobachevsky 
plane. The key result of our study consists in the temporary gauge invariance of 
the billiard representation, once provided the use of very general Misner-Chitre- 
like variables. 



1 Introduction 



Since Belinski-Kalatnikov-Lifshitz (BKL) derived the oscillatory regime which 
characterizes the behaviour of the Bianchi type VIII and IX cosmological models 
[0 - (the so-called Mixmaster universe ||) near a physical singularity, a wide 
literature faced over the years this subject in order to provide the best possible 
understanding of the resulting chaotic dynamics. 

The research activity developed overall in two different, but related, directions: 
on one hand the dynamical analysis was devoted to remove the limits of the 



BKL approach due to its discrete nature (by analytical treatments Ig-JT^l and 



numerical simulations |TTJ-||14|), on the other one to get a better characterization 



of the Mixmaster chaos (especially in view of its properties of covariance [15|- [18| ) 



The first line of investigation provided satisfactory representations of the Mix- 
master dynamics in terms of continuous variables (leading to the construction of 
an invariant measure for the system ||, ||), as well as detailed numerical descrip- 
tions (allowing to make precise validity tests on the obtained analytical results) 
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The efforts to develop a precise characterization of the chaoticity observed in 
the Mixmaster dynamics found non-trivial difficulties due to the impossibility, 
or in the best cases the ambiguity, to apply the standard chaos indicators to 
relativistic systems. In particular, the puzzle consisting of numerical simulations 
which were providing, zero |T^] and non-zero ]T7|] Lyapunov exponents by using 
different time variables, has been solved by realizing the non-covariant nature 
of these indicators and their inapplicability to hyperbolic manifolds The 
existence of these difficulties prevented, up to now, to say a definitive word about 
the covariance of the Mixmaster chaos, with particular reference to the possibility 
of removing the observed chaotic features by a suitable choice of the time variable 
(apart from the indication provided by 0). 

Indeed a valuable framework of analysis of the Mixmaster evolution, able to 
join together the two considered points of view, relies on an Hamiltonian treat- 
ment of the dynamics in terms of Misner-Chitre-like variables 12IJ . This formula- 
tion allows to individualize the existence of an asymptotic (energy-like) constant 
of motion when performed an ADM reduction. By this result the stochasticity 
of the Mixmaster can be treated either in terms of the statistical mechanics (by 
the microcanonical ensemble), either by its characterization as isomorphic to a 
billiard on the two-dimensional Lobachevsky space \TL\. 



The aim of this work is to show how the representation of the Mixmaster 
dynamics as a "stochastic scattering" isomorphic to a billiard on the Lobachevsky 
space can be constructed independently of the choice of a time variable, simply 
providing very general Misner-Chitre-like coordinates. Up to the limit of the 
adopted approximation on the form of the potential term, our analysis shows, 
without any ambiguity, that the Mixmaster stochasticity can not be removed by 
any redefinition of the time variable. 
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More precisely in Section 2 we construct the Hamiltonian formulation of the 
Mixmaster dynamics, which in Section 3, by a standard Arnowitt-Deser-Misner 
(ADM) reduction and the study of the asymptotic form of the potential, allows 
to individualize an appropriate energy-like constant of motion (nothing more 
then the ADM kinetic energy) in any temporal gauge. In Section 4 the reduced 
Hamiltonian principle is rewritten in terms of a geodesic one on the induced 
Jacobi metric. Finally in Section 5 we derive the statistical implications due to 
the properties of geodesic flow and provide brief concluding remarks. 



2 The Hamiltonian Formulation 

The geometrical structure of the Bianchi type VIII and IX spacetimes, i.e. of the 
so-called Mixmaster Universe models, is summarized by the line element [ffj 

ds 2 = —N(i]) 2 dr] 2 + e 2a (e 213 ),^ a'a? (1) 

where N(rj) denotes the lapse function, a 1 are the dual 1-forms associated with the 
anholonomic basis [] and /3y is a traceless 3x3 symmetric matrix diag(/3n, /? 22 , ^33); 
a, N, flij are functions of r\ only. Parameterizing the matrix /3y by the usual Mis- 
ner variables 0] 

fti = (3 + + V3(3. 

f3 22 = P+-V3/3- (2) 
/9 33 = ~2f3 + 

the dynamics of the Mixmaster model is described by a canonical variational 
principle 

5I = djLdr] = J (3) 

with Lagrangian L 



L- 6D 
~ AT 



-a' 2 + /V 2 + P-' 2 \ - (a, p+, /?_) . (4) 



N, 



1 The dual 1-forms of the considered models are given by: 



(Bianchi VIII): 



(Bianchi IX): 



a 1 = 


— sinh tj) sinh 6d(f> 


+ cosh tpdO 


a 2 = 


— cosh ip sinh 9d(f> 


+ sinh tpdO 


a 3 = 


cosh9d(f> + dijj 




a 1 = 


sin ip sin 9d<j> + 


cos 4>d9 


a 2 = 


— cos ijj sin Odcj) - 


- sin-0ii0 


a 3 = 


cos 6d(j) + dip 
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Here ()' = -f, D = det e a+ft -? = e 3a and the potential V (a, /3+, /?_) reads 



V =- ^D 4Hl + D 4H2 + D 4Ha ^j - D 2Hl+2H2 ± D 2H ' 2+2Hs ± D 2Hs+2Hl , (5) 

where (+) and (— ) refers respectively to Bianchi type VIII and IX, and the 
anisotropy parameters Hi (i = 1,2,3) denote the functions || 



Hi 

H 2 
H 3 



1 | (3 + + V3p. 
3 3a 
1 | P+-V3/3- 
3 3a 

3 3a 



(6) 



In the limit D — ► the second three terms of the above potential turn out to be 
negligible with respect to the first one. Let's introduce the new (Misner-Chitre- 
like) variables 



a 



_ e /(r) C 



J(r) 



J(r) 



£ 2 - 1 COS 9 



(7) 



£ 2 - 1 sin 9 . 



with / denoting a generic functional form of r, 1 < £ < oo and < 9 < 2ir. 
Then the Lagrangian (|J) reads 



QD 
It 



M + (^)"( 4 ._ 1 )_ (e0 



l2 



N 

~ ^V(/(r),e, 



In terms of / (r), £ and 6* we have 

D = exp{-3£-e /(r) } 



(9) 



and since D — > toward the singularity, independently of its particular form, in 
this limit / must approach infinity The Lagrangian @ leads to the conjugate 
momenta 



Pr 



Pi 



dL 




d? ~ 


N { 


dL 


12D e 2 f 


d?~ 


N £ 2 - 1 


dL 
d9' ~ 





dr 



(10) 
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which by a Legendre transformation make the variational principle (|3]) assume 
the Hamiltonian form 



r ( Ne~ 2f \ 

s J Uc + P9 0' + py - -^-nj d v = o 



(11) 



being 



n = + (i 2 - l) + + 24 W . (12) 



dr 



3 Reduced Variational Principle 



By variating (|TT|) with respect to iV we get the constraint TC — 0, which solved 
provides 

df „ df 



-p T = -jL. H ADM = ■ \]e 2 + 2AVe 2 f (13) 
where 

in terms of which the variational principle (|ll|) reduces to 

S J (p£ + PeO' - f"H A DM) d V = 0. (15) 
Since the equation for the temporal gauge actually reads 

N (V) = #^fr' , (16) 

it ADM OT 

our analysis remains fully independent of the choice of the time variable until the 
form of / and r' is not fixed. 

The variational principle ( |T5| ) provides the Hamiltonian equations for £' and 

e = ^(e-i)p, 

it ADM 

9 ' = ^wle^ry (17) 

Furthermore can be straightforward derived the important relation 

d (Hadm/') d {HadmI') d (HadmD 9 (Hadm/') n „x 



dr/ dr] df df 



2 In this study the corresponding equations for p'^ and p' g are not relevant. 
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i.e. explicitly 

In this reduced Hamiltonian formulation, the functional / (77) plays simply the 
role of a parametric function of time and actually the anisotropy parameters Hi 
(i = 1, 2, 3) are functions of the variables £, 9 only 



Hi = \- 1 (cos 9 + VSsmO 

1 /£2 ~T" 

#2 = 3- 3^ (cos vising) (20) 

# 3 = i + 2 — ~ - cosfl. 
3 3 3£ 

Finally, toward the singularity (D — > i.e. / — > 00) by the expressions ([5], ^ 
, we see that f] 

d Z = (eM • (21) 



5/ 

Since in the domain r#, where all the Hi are simultaneously greater than 0, the 
potential term U = e 2 *V can be modeled by the potential walls 



U 00 = Goo (#!(£, 0)) + 600(^2(^0)) + 600(^3^,^)) (22) 

J +00 £/ x < 
I i/ x > 



6oo fx) 



therefore in Th the ADM Hamiltonian becomes (asymptotically) an integral of 
motion 



WJt fl , c r / ^adm = v^ 2 + 2A-U = e = E = const. 
V{£, fl}er H i 8H AnM _ dE _ n (23) 
la/ a/ u - 

The key point for the use of the Misner-Chitre-like variables relies on the 
independence of the time variable for the anisotropy parameters Hi. 



4 The Jacobi Metric and the Billiard Represen- 
tation 

Since above we have shown that asymptotically to the singularity (/ — > 00, i.e. 
a — > —00) dTiADNi/df = i.e. TCadm = € = E = const., the variational principle 

3 By 0() we mean terms of the same order of the inclosed ones. 
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( |15|) reduces to 

^ / (pe<% + Pedd - Edf) = 6 J (p € d£ + Pe d9) = , (24) 

where we dropped the third term in the left hand side since it behaves as an exact 
differential. 

By following the standard Jacobi procedure |^TJ to reduce our variational 
principle to a geodesic one, we set x al = g ab Pb, and by the Hamiltonian equation 
( |i~7l) we obtain the metric 

By these and by the fundamental constraint relation 

:e-i)p, 2 +f^ = E\ (26) 



g ab x a 'x b ' = f -Ue- lW + #t) = f' E - (27) 



we get 

By the definition x a/ = = w a J, the (0) rewrites 

drj 

which leads to the key relation 



9atu"u b (^\ =f'E, (28) 



*n = f-^is. (29) 

Indeed the expression ( p9|) together with + = allows us to put the 
variational principle in the geodesic form: 



5 J f'Ed V = 5 J ^g ab u-u b f'E ds = 6 J ^G ab u-u b ds = (30) 

where the metric G a b = f'Eg a b satisfies the normalization condition G a bU a u b = 1 
and therefore Q 

ds „, ds , . 

= = (31) 



4 We take the positive root since we require that the curvilinear coordinate s increases mono- 
tonically with increasing value of /, i.e. approaching the initial cosmological singularity. 
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Summarizing, in the region IV the considered dynamical problem reduces to a 
geodesic flow on a two dimensional Riemannian manifold described by the line 
element 

f 2 _ 1 + ^ 

Now it is easy to check that the above metric has negative curvature, since the 



ds 2 = E 



Dde 2 



(32) 



associated curvature scalar reads R 



E 2 ' 



therefore the point-universe moves 



over a negatively curved bidimensional space on which the potential wall (|5|) 
cuts the region Th- By a way completely independent of the temporal gauge we 
provided a satisfactory representation of the system as isomorphic to a billiard 



on a Lobachevsky plane 21 



5 Invariant Lyapunov Exponent 

In order to characterize the dynamical instability of the billiard in terms of an 
invariant treatment (with respect to the choice of the coordinates £, 9), let us 
introduce the following (orthonormal) tetradic basis 

w< = (°-17F=t) ' (33) 

Indeed the vector v l is nothing else than the geodesic field, i.e. 

DV ' dv ' + r kl vV = , (34) 



ds ds 

while the vector w l is parallely transported along the geodesic, according to the 
equation 

^ = f^ + r>V = o, (35) 

ds ds kt K ' 

where by T l kl we denoted the Christoffel symbols constructed by the metric (53). 
Projecting the geodesic deviation equation along the vector w l (its component 
along the geodesic field v 1 does not provide any physical information about the 
system instability), the corresponding connecting vector (tetradic) component Z 
satisfies the following equivalent equation 

d 2 Z Z 

7?-W (36) 

This expression, as a projection on the tetradic basis, is a scalar one and therefore 
completely independent of the choice of the variables. Its general solution reads 

Z (s) = c\e® + C2e~E , c± t 2 = const. , (37) 
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and the invariant Lyapunov exponent defined as p2 



= sup lim — , (38) 



in terms of the form (|3"7|) takes the value 



A„ = i>0. (39) 

When the point-universe bounces against the potential walls, it is reflected from 
a geodesic to another one thus making each of them unstable. Though up to the 
limit of our potential wall approximation, this result shows without any ambiguity 
that, independently of the choice of the temporal gauge, the Mixmaster dynamics 
is isomorphic to a well described chaotic system. Equivalently, in terms of the 
BKL representation, the free geodesic motion corresponds to the evolution during 
a Kasner epoch and the bounces against the potential walls to the transition 
between two of them. By itself, the positive Lyapunov number (|39D is not enough 
to ensure the system chaoticity, since its derivation remains valid for any Bianchi 
type model; the crucial point is that for the Mixmaster (type VIII and IX) the 
potential walls reduce the configuration space to a compact region (T H ), ensuring 
that the positivity of A„ implies a real chaotic behaviour (i.e. the geodesic motion 
fills the entire configuration space). 

Summarizing, our analysis shows that for any choice of the time variable, 
we are able to give the above stochastic representation of the Mixmaster model, 
provided the factorized coordinate transformation in the configuration space 

a =-e /{r) a(#,0 

13+ = e^b + (6,0 (40) 
/?_ = e'< r >&_(0,£) , 

where /, a, b+ denote generic functional forms of the variables r, 9, £. 

It is worth noting that the success of our analysis, in showing the time gauge 
independence of the Mixmaster chaos, relies on the use of a standard ADM re- 
duction of the variational principle (which reduces the system by one degree of 
freedom) and overall because, adopting Misner-Chitre-like variables, the asymp- 
totic potential walls are fixed in time. The difference between our approach and 
the one presented in |23| (see also for a critical analysis |24|) consists effectively in 



these features, though in those works is even faced the problem of the Mixmaster 
chaos covariance with respect to the choice of generic configuration variables. 

We are very grateful to Remo Ruflini for his valuable comments on this sub- 
ject. 
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